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The recently discovered (Lio.sFeo.2 )OHFeSe superconductor provides a new platform for exploiting 
the microscopic mechanisms of high-Tc superconductivity in FeSe-derived systems. Using density 
functional theory calculations, we first show that substitution of Li by Fe not only significantly 
strengthens the attraction between the (Lio. 8 Feo. 2 )OH spacing layers and the FeSe superconducting 
layers along the c axis, but also minimizes the lattice mismatch between the two in the ab plane, 
both favorable for stabilizing the overall structure. Next we explore the electron injection into 
FeSe from the spacing layers, and unambiguously identify the Feo .2 components to be the dominant 
atomic origin of the dramatically enhanced interlayer charge transfer. We further reveal that the 
system strongly favors collinear antiferromagnetic ordering in the FeSe layers, but the spacing layers 
can be either antiferromagnetic or ferromagnetic depending on the Feo .2 spatial distribution. Based 
on these understandings, we also predict (Lio. 8 Coo. 2 )OHFeSe to be structurally stable with even 
larger electron injection and potentially higher Tc. 

PACS numbers: 74.70.-b, 74.62.Dh, 74.25.Ha, 84.30.Bv 


Recently, two groups discovered (Lii_a,Fea;)OHFeSe 
{x ^0.2) as a new class of superconductors with high su¬ 
perconducting transition temperatures (Tc), and demon¬ 
strated coexistence of superconductivity and antiferro¬ 
magnetism (AFM) or ferromagnetism (FM) in these sys¬ 
tems [1, 2]. Such systems provide several appealing fea¬ 
tures for potentially revealing the likely superconducting 
mechanisms [1-4]. However, despite quite a few subse¬ 
quent experimental studies, limited knowledge has hith¬ 
erto been obtained about the likely dominant roles played 
by the spacing layers [3-7], especially on how essential the 
dopant Fex atoms are in establishing superconductivity. 
Intriguing questions include the magnetic ordering of the 
spacers, which appears to be inconsistent between the 
two pioneering experiments [1, 2]. More importantly, it 
is critically desirable to understand the structural sta¬ 
bility of the systems formed by the weakly interacting 
spacers and superconducting FeSe monolayers [1, 2], as 
well as the atomic origin of the electron doping from the 
(Lii_a,Fea^)OH layers. Insights into those important is¬ 
sues, in turn, will be instrumental in searching for new 
FeSe-based superconductors with potentially higher Tc. 

In this Letter, we use density functional theory (DFT) 
calculations to investigate the dominant roles of the 
spacers in establishing the high-Tc superconductivity 
of (Lio. 8 Feo. 2 )OHFeSe, with particular emphasis on the 
Feo .2 atoms. We find that substitution of Li by Fe 
strengthens the structural stability both in the ab plane 
and along the c axis. By further exploring the charge 
transfer, we identify the Feo .2 atoms to be the atomic ori¬ 
gin of significant electron injection into FeSe. In addition, 


we obtain the ground-state magnetic order, and explain 
the seemingly controversial experimental observations on 
the magnetic ordering of the spacers. Based on these 
findings, we also predict a stable (Lio. 8 Coo. 2 )OHFeSe 
structure with larger electron doping into FeSe, poten¬ 
tially resulting in an even higher Tc- 



FIG. 1. (a) A side view of the 1x1 (Li/Fe)OHFeSe, and 
top views of (b) the FeSe and (c) (Lio. 8 Feo. 2 )OH layer in a 
\/T0 X \/Td cell. The green shaded area is the 1x1 cell, while 
the yellow is \/5 x which is compatible with states of 
nonmagnetic (NM), FM, or checkerboard AFM in FeSe. The 
collinear AFM in FeSe requires an expanded cell of a/TO- 

We first perform lattice-constant relaxation of 
(Lio. 8 Feo. 2 )OHFeSe with different functionals and van der 
Waals (vdW) approaches, in order to choose the proper 
theoretical framework that can describe the current sys¬ 
tems [8]. In constructing the supercells, we have two 
considerations: ( 1 ) (Lio.8Feo.2)OH layers have an occu¬ 
pation of Li/Fe with a ratio of ^4 [1], and thus the small¬ 
est cell is X ^/ 5 ; (2) as discussed explicitly below, 
the ground-state magnetism is obtained to be collinear 
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AFM for FeSe, consistent with several previous stud¬ 
ies of FeSe/STO [9-11]. Therefore, a x ^/l0 cell 
is used to optimize the lattice parameters of ground- 
state (Lio. 8 Feo. 2 )OHFeSe. Figs, lb and Ic display the 
atomic structures of the FeSe and (Lio. 8 Feo. 2 )OH layers; 
in a first study, the Feo .2 atoms in the spacer are or¬ 
derly distributed and thus form a square lattice, where 
the nearest-neighboring Fe-Fe distance is a constant. Our 
detailed studies show that the clustering or disordering 
of the Feo .2 atoms is energetically less favorable [8]. 

TABLE 1. Calculated lattice constants of (Lio. 8 Feo. 2 )OHFeSe 
and LiOHFeSe. The vdW methods are based on the PBE 
functional. Numbers are in units of A. These results are ob¬ 
tained from optimization of the \/Tb x \/l0 cells, and the data 
is normalized to the 1x1 cells. 


(Lio. 8 Feo. 2 )OHFeSe 

LiOHFeSe 

EXP^ 

non-vdW 

vdW 

non-vdW 

LDA 

PBE 

DFT-D2 

DFT-TS 

vdW-DF2 

PBE 

a = 3.786 

FTE^ 

3.79 

3.72 

3.69 

FTE‘> 

3.68 

c = 9.288 

< 8.20 

9.38 

8.52 

8.30 

> 9.78 

11.16 


^ EXP: experimental data at room temperature [1]. 
b ptE: fix to the experimental parameter. 


The optimized parameters are shown in Table I, cal¬ 
culated from the energy dependence on the lattice con¬ 
stants. For LDA and vdW-DF2, we freeze a to the ex¬ 
perimental data [1], and find a large deviation in c. We 
therefore relax only c even without convergence for these 
two columns, while for the rest, both a and c are op¬ 
timized until convergence is well reached. Between the 
non-vdW approaches, LDA significantly underestimates 
the lattice constant, while PBE slightly overestimates the 
axes, as commonly observed for the two functionals. The 
implemented vdW corrections as the empirical pairwise 
forms of Cq/Rq [12, 13] in DFT-D2 and DFT-TS are 
found to over-reduce c; more surprisingly, the non-local 
vdW functional [14-16] included in vdW-DF2 is found 
to severely overestimates c. Here the failures of these 
vdW methods may originate from the excessively large 
Cq coefficients due to the neglect of screening effect and 
many-body formalism for the former two, and from the 
imprecise exchange functional characterized in the lat¬ 
ter. These problems call for in-depth investigations and 
possibly require major improvements on the methodol¬ 
ogy aspect. Presently, given the relatively small errors 
{a: +0.16%, c: +0.96%), we conclude that the PBE 
functional without including either of the popular vdW 
corrections yields a reasonable description, and thus the 
following results are all based within this framework. 

To determine the roles of Feo .2 in the structural proper¬ 
ties, we calculate the lattice constants of LiOHFeSe. The 
Fe atoms are removed from the spacer, and the collinear 
AFM order is found to remain for the FeSe layer. The 
obtained parameters (Table I) are reduced by 3.01% in a 
and enlarged by 19.01% in c individually, from the lat¬ 
tice constants of (Lio. 8 Feo. 2 )OHFeSe calculated by PBE. 


In the following, we analyze the effects on a and c sepa¬ 
rately in more details. 

Experimentally, FeSe is inferred to be compressed in 
the ab plane of (Lio. 8 Feo. 2 )OHFeSe [1], where intuitively 
the strain should be caused by the relatively smaller lat¬ 
tice of (Lio. 8 Feo. 2 )OH. Based on our results, the Feo .2 
atoms indeed have expanded the ab lattice; otherwise, 
the spacer could be even smaller in ab, and the corre¬ 
sponding even larger mismatch with FeSe would make 
it unlikely to form a stable intercalated structure. To 
confirm this conjecture, we further calculate the lattice 
parameters of FeSe, (Lio. 8 Feo. 2 )OH, and LiOH monolay¬ 
ers separately, whose a values are found to be 3.74 A, 
3.73 A, and 3.59 A, respectively. These values, together 
with that of bulk (Lio. 8 Feo. 2 )OHFeSe, suggest that the in¬ 
plane compressed nature [1] of FeSe may be inaccurate; 
however, our conjecture remains valid, given the small {or 
large) lattice mismatch between FeSe and (Lio. 8 Feo. 2 )OH 
{or LiOH). Therefore, the Feo .2 atoms play a vital role in 
minimizing the lattice mismatch between the spacer and 
FeSe, which enables the formation of the commensurate 
stacking as shown in Fig. la. 

The FeSe and (Lio_sFeo_ 2 )OH layers are both in-plane 
stretched slightly to 3.79 A when they are stacked al- 
ternatingly. This seems counter-intuitive, as the con¬ 
tacting interface generally adopts an intermediate lat¬ 
tice constant. To understand this observation, we ex¬ 
amine the structures by calculating the thickness d of 
each layer, defined as the distance from the upper Se 
or H to the lower Se or H position along the c direc¬ 
tion for the FeSe monolayer or (Lio. 8 Feo. 2 )OH spacing 
layer. We find d is quenched from 3.48 to 3.33 A for 
(Lio. 8 Feo. 2 )OH, while expanded from 2.89 to 2.91 A for 
FeSe, when the two layers are alternatively assembled 
into bulk. The volume expansion (from 3.74^x2.89 to 
3.79^x2.91 A^) of FeSe suggests an increased Coulomb 
repulsion internally, which strongly indicates electron in¬ 
jection. The (Lio. 8 Feo. 2 )OH layer is found to be rippled, 
with an amplitude of ^0.12 A for H deviating from their 
average position along the c-axis, making its in-plane lat¬ 
tice relatively easy to be stretched by interacting with the 
expanded FeSe, and the thickness tends to decrease as a 
compensation. Combining all the dimensions, the volume 
shrinks (from 3.73^x3.48 to 3.79^x3.33 A^), supposed to 
be caused by electron depletion. To explain the counter¬ 
intuitive observation, here we emphasize the subtle cor¬ 
relation between lattice expansion and electron doping, 
which probably has been neglected in other FeSe-based 
systems (such as FeSe/STO and alkali-intercalated FeSe) 
that rely on charge transfer to enhance Tc. 

We now examine specifically how Feo .2 influences the 
structures along the vertical directions. The c axis is re¬ 
duced from ll.16Ato9.38A when the Feo .2 atoms are in¬ 
corporated into the spacer. Such a dramatic change sug¬ 
gests a greatly enhanced attraction between the spacer 
and FeSe. We thus calculate the cohesive energy: Ec = 
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EpeSe^^Spacer -Esuik, where EpeSe IS the total energy 
of the FeSe monolayer, Espacer is that of (Lio. 8 Feo. 2 )OH 
{or LiOH), and Esuik is that of the combined bulk. 
EpeSe and Espacer are calculated using the in-plane lat¬ 
tice of the corresponding bulk system. Our results show 
that, Ec per formula unit is 0.02 eV for LiOHFeSe and 
0.35 eV for (Lio. 8 Feo. 2 )OHFeSe. For the former, the value 
is even smaller than Ec per carbon in graphite [17], and 
the average bond distance d of the nearest Se-H is 3.50 
A, clearly in the regime of weak vdW interaction. In 
contrast, for (Lio. 8 Feo. 2 )OHFeSe, Ec is stronger than 
that of vdW interaction, but still weaker than the typical 
strength of chemical bonding. This moderate interaction 
originated from the incorporated Feo .2 substantially sta¬ 
bilizes the structure vertically, while still allows mechan¬ 
ical cleavage of the crystal [4, 6, 7]. More importantly, 
such an increased Ec and the closer contact {d = 3.11 
A) between FeSe and (Lio. 8 Feo. 2 )OH collectively enable 
enhanced charge transfer between the layers. 

Following the above structural indications, we now in¬ 
vestigate the detailed nature of charge transfer. We cal¬ 
culate the charge density difference Ap between the com¬ 
bined (Lio. 8 Feo. 2 )OHFeSe bulk system and the sum of the 
isolated FeSe and (Lio. 8 Feo. 2 )OH layers. To have a quan¬ 
titative picture, we plot the plane-averaged Ap along the 
c axis (Ap;^) in Fig. 2a. To understand the roles of Feo. 2 , 
we also calculate /S.pz for LiOHFeSe at the fixed lattice 
constant of (Lio. 8 Feo. 2 )OHFeSe and at the relaxed lattice, 
as shown in Figs. 2b and 2c respectively. 
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FIG. 2. Ap^ of (a) (Lio. 8 Feo. 2 )OHFeSe and (b,c) LiOHFeSe at 
the lattice constant of (Lio. 8 Feo. 2 )OHFeSe and of LiOHFeSe, 
(d) (Lio. 8 Mno. 2 )OHFeSe, and (e) (Lio. 8 Coo. 2 )OHFeSe at their 
own relaxed lattices. These results are calculated from the 
ground states of \/l0 x \/l0 cells, and have been normalized 
to 1 X 1 cells. The red and blue shaded areas indicate electron 
accumulation and depletion, respectively. 

We find that, when the (Lio. 8 Feo. 2 )OH and FeSe layers 
merge in bulk, the charge density between Fe and Se in 
FeSe has an electron-depletion region, and this amount 
of charge is mostly transferred to the Se planes to form 


interactions with the spacers. Such features are not ob¬ 
served in LiOHFeSe, where the interlayer Ec is much 
weaker. In each plot, we shift the average 2 ; position of 
the Fe atoms in FeSe to be located exactly at the cen¬ 
ter of the supercell. The charge redistribution is mirror 
symmetric with respect to the middle line in Figs. 2a 
and 2b. However, in Fig. 2c, the layers are loosely stacked 
with the interlayer distance found to be slightly alternat¬ 
ing (the distance between the Se and H planes differs by 
^0.15 A at the two opposite sides of FeSe), making the 
curve not to be precisely symmetric. To see the dop¬ 
ing levels of FeSe, we integrate /S.pz over 2 ; within the Se 
boundary on each side. We define a boundary between H 
and Se, so that the distance from Se (H) to the boundary 
is proportional to the atomic radius of Se (H) [18]. The 
electron injection p/ into FeSe is calculated to be 0.051 
electrons per FeSe in (Lio. 8 Feo. 2 )OHFeSe. For LiOHFeSe, 
pi = 0.015e/Fe and 0.005e/Fe for structures with the 
(Lio. 8 Feo. 2 )OHFeSe and LiOHFeSe lattice constant, re¬ 
spectively. The estimated value of 0.051 is close to the 
experimentally measured value of ^0.08 e/Fe [3], and 
the difference might be due to the absence of structural 
defects or disordering of Feo .2 in the supercell systems. 

From a broader perspective, a sufficiently large elec¬ 
tron doping is essential to realize high Tc in FeSe-based 
superconductors. In the FeSe/STO systems, O vacancies 
in STO are believed to be the atomic origin of large in¬ 
terlayer charge transfer [19]. Here in (Lio. 8 Feo. 2 )OHFeSe, 
one may naturally think that the atomic origin of elec¬ 
tron injection is the Feo .2 atoms in the spacer, because 
of their generally higher oxidation state than Li. Indeed, 
our results agree with this conjecture, and demonstrate 
that the Feo .2 atoms enhance the electron injection into 
FeSe in two ways. First, by comparing Figs. 2a and 2b 
where the lattice constants are identical, we clearly see 
the large contribution of the Feo .2 d orbitals to the charge 
transfer when Li is substituted. Second, the increased 
doping in Fig. 2b compared to Fig. 2c indicates that the 
closer interlayer coupling caused by Feo .2 also boosts the 
amount of charge transfer. These two effects are both 
crucial, and should be instrumental in search for other 
novel spacer-intercalated superconductors. 

The interlayer charge transfer results in the higher Tc 
in bulk (Lio. 8 Feo. 2 )OHFeSe than in bulk FeSe; however, 
the doping level of FeSe in (Lio. 8 Feo. 2 )OHFeSe is still 
lower than that in FeSe/STO [3], suggesting the feasi¬ 
bility of other structural design to further enhance Tc. 
Based on the above studies of the roles of Feo. 2 , we 
can attempt to substitute Li by other elements X in¬ 
stead of Fe in the spacers, to see if such elements can 
also stabilize the structure, and more importantly, in¬ 
duce larger charge transfer. Presently, we have examined 
X = Mn or Co, which potentially has higher oxidation 
state than that of Fe. Our calculations show that both 
(Lio. 8 Mno. 2 )OHFeSe and (Lio. 8 Coo. 2 )OHFeSe are indeed 
structurally stable, with their lattice constants similar to 
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that of (Lio. 8 Feo. 2 )OHFeSe. Furthermore, based on the 
charge transfer shown in Figs. 2d and 2e, pi is calculated 
to be 0.044e/Fe for (Lio. 8 Mno. 2 )OHFeSe and 0.060e/Fe 
for (Lio. 8 Coo. 2 )OHFeSe. The Mn and Co atoms do not 
appear to contribute to the charge injection of FeSe as 
much as Fe, but significantly facilitate more electron 
transfer from the O atoms in the spacer. The depleted 
area in FeSe of (Lio. 8 Coo. 2 )OHFeSe is also found to be 
smaller than that in (Lio. 8 Feo. 2 )OHFeSe. The larger pi 
of (Lio. 8 Coo. 2 )OHFeSe suggests a higher Tc than that of 
(Lio. 8 Feo. 2 )OHFeSe. More studies should be made along 
this line of structural design, involving other possible el¬ 
ements and substitution concentrations > 0.2, to further 
identify the best candidates in this family of materials to 
realize much higher 

We next focus on the magnetic properties, as mag¬ 
netism is generally related to superconductivity in the 
FeSe-derived systems. We examine four different mag¬ 
netic orders in each of the FeSe and (Lio. 8 Feo. 2 )OH lay¬ 
ers, including NM, FM, checkerboard AFM, and collinear 
AFM, and investigate their possible combinations for in¬ 
terlayer coupling. Our calculations show that FM is 
unable to be established in FeSe, as well as NM in 
(Lio. 8 Feo. 2 )OH. The computed magnetic moment of Feo .2 
is about 3.53/iB, and for that of Fe in FeSe, it is ^ 2.25/iB 
in the collinear AFM and ^ 1.86/iB in the checkerboard 
AFM state. The larger moment of Feo .2 suggests that 
the spin magnitude is reduced in a closer packing of the 
Fe atoms in FeSe. By comparing the total energies of 
the different orders, we obtain the magnetic ground state 
of (Lio_sFeo, 2 )OHFeSe (Fig. 3). Both layers exhibit a 
collinear AFM order in the Fe square lattices, and the 
interlayer spins are aligned parallel to each other. AFM 
is calculated to be slightly more stable than FM by only 
2 meV per Fe in (Lio. 8 Feo. 2 )OH. Such results of the frag¬ 
ile AFM ground state in the spacer also help to explain 
the recent experimental observations [20]. In addition, 
for the newly designed systems, (Lio. 8 Coo. 2 )OHFeSe has 
the same spin configuration as (Lio. 8 Feo. 2 )OHFeSe, while 
(Lio. 8 Mno. 2 )OHFeSe exhibits anti-parallel alignment be¬ 
tween the collinear AFM FeSe and (Lio. 8 Mno. 2 )OH lay¬ 
ers. 



FIG. 3. (a) Ground-state magnetism of (Lio. 8 Feo. 2 )OHFeSe 
in the (b) FeSe and (c) (Lio. 8 Feo. 2 )OH layers. Elements other 
than Fe have moments of ~0, and are colored in gray. In 
(a), the arrows spanning over the Fe rows indicate the same 
spin orientation in the row. The green lines in (b) and (c) 
connect the nearest or next-nearest Fe neighbors, whose mag¬ 
netic couplings are considered in the Heisenberg model. 


Based on these DFT inputs, we quantitatively analyze 
the coupling strengths in (Lio. 8 Feo. 2 )OHFeSe using the 
Heisenberg model on square lattices [11, 21, 22]. We use 
an approximate Hamiltonian: 

<ij> <^ij^ <ij> <^ij^ <ij> 

( 1 ) 

The three terms respectively account for couplings in 
FeSe, in (Lio. 8 Feo. 2 )OH, and between the two layers. 
<ij> and <^ij^ represent summations over the near¬ 
est and next-nearest neighbors. Based on the energies 
of different magnetic states, we estimate the coupling 
strengths (Table II) [8]. We note that the computed Jn 
and Ji 2 in FeSe are quite close to previous results [11, 21]. 

TABLE IF Galculated J values of (Lio. 8 Eeo. 2 )OHEeSe, in 
which the Eeo .2 atoms are uniformly distributed in spacer. 



Jii 

J 12 

J 21 

J 22 

J 3 

Unit: meV 

83.18 

47.82 

0.05 

-0.02 

-1.94 


Generally speaking, the order of the coupling strength 
is: ^ 11,12 ^ Js ^ '^ 21,22 5 just the reverse of the Fe- 
Fe pairing distance (diiq 2 < ds < ^ 21 , 22 )- The cou¬ 
plings in (Lio. 8 Feo. 2 )OH are quite weak, making its mag¬ 
netism more dominated by the relatively larger coupling 
between the layers. The magnetism of FeSe and the neg¬ 
ative J 3 thus result in a collinear AFM ground state in 
(Lio. 8 Feo. 2 )OH (Fig. 3c), with the spin direction of each 
site parallel to that of the corresponding Fe in FeSe. Fur¬ 
thermore, the much weaker J 3 compared to Jn and J 12 
suggests that, Feo .2 should play a minimal role in directly 
influencing the magnetism of FeSe by magnetic interplay 
between the layers. This fact also indicates that the mag¬ 
netic coupling between FeSe and the spacer is unlikely to 
play an important role, or probably even undesirable in 
establishing high Tc. 

It is notable that, the above ground-state magnetism 
is computed using an ordered distribution of Feo .2 iu the 
spacer. In the actual systems, certain degrees of disor¬ 
dering is unavoidable. This structural fact could pos¬ 
sibly change the coupling strengths in the spacer and 
the Fe sites in FeSe that are coupled closely to the Feo .2 
atoms. By calculating the ground-state magnetism, here 
we demonstrate that disordering of Feo .2 can possibly re¬ 
sult in a FM order in the spacer [8] . This finding may help 
to clarify the controversial experimental observations of 
different magnetic orders in (Lio.8Feo.2)OH [ 1 , 2 ]. 

Overall, despite several limitations [8], this study has 
revealed the dual roles of Feo .2 in the structural stabil¬ 
ity and electronic charge injection into FeSe. Such roles 
are critically important in the fabrication of FeSe-based 
high-Tc superconductors, and may provide new insights 
into exploration of the likely pairing mechanisms. Our 
predicted (Lio. 8 Coo. 2 )OH superconductors with larger 
charge transfer and potentially higher Tc also calls for 
experimental fabrication and validation. 
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